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Abstract 

Substrate disorder effects on the scaling properties of growing crystalline surfaces in solidification 
or epitaxial deposition processes are investigated. Within the harmonic approach there is a phase 
transition into a low-temperature (low-noise) superrough phase with a continuously varying dynamic 
exponent z > 2 and a non-linear response. In the presence of the KPZ nonlinearity the disorder causes 
the lattice effects to decay on large scales with an intermediate crossover behavior. The mobility of the 
rough surface has a complex dependence on the temperature and the other physical parameters. 
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Much progress has been achieved recently in the understanding of kinetic roughening in 
nonequihbrium surface growth. In particular spatial and temporal scaling properties have 
been predicted for surfaces of crystals growing either by solidification or from epitaxial deposition.^ 
The simplest analytical models neglect the discrete structure of the crystal and the scaling prop- 
erties are determined from the university classes of the continuous kinetic equations which govern 
the growth process. 

The lattice structure^"^^ has, however, crucial effects on the behavior with a kinetic phase 
transition mirroring the equilibrium roughening transition (ERT). The transition is governed 
by the noise, i.e. thermal and from the inherent stochasticity, of the growth process. In 
presence of a strong noise the surface is rough (in analogy with the rough equilibrium phase for 
T > Tr, where is the ERT temperature). In this phase the lattice effects are unimportant and 
the dynamic properties are the same as in the respective continuum models. In the low-noise 
phase the surface is smooth on an intermediate scale (which becomes larger the smaller the 
average growth rate is) with drastically different dynamic properties. ^'^ 

In view of the crucial importance of the discrete structure at the low-noise regime, the fol- 
lowing question should be addressed: What if the substrate on which the solid grows is not 
perfectly smooth? While all previous studies have assumed a perfectly flat substrate, in many 
potential realizations that will not be the case. The effects of quenched disorder in the substrate 
on the kinetic scaling properties of the growing crystalline surface are the subject of the present 
letter. As detailed below we find that in presence of substrate disorder the low-noise (or low 
temperature) regime has dramatically different kinetic properties. 

The scaling properties are manifested in the height-height correlation function: 

C{L,t) = {Ah'{L,T)) = {[h{x + L,t + T) - h{x,t)]') ^ L'y{T/L^), (1) 
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h{x, t) is the height of the surface at point x and time t. The roughness exponent a characterizes 
the long-time hmit of the self-affine fractal structure as represented by the "surface width" 
w{L,t) = (A/i^(L, r))^^^ fti L". At early stage of the growth (r <^ L^) w where /? = a/z 
and z is the dynamic exponent. Theoretical studies of kinetic surfaces have been conducted 
within two approaches: the "harmonic" approach^'^'''' and the "non-linear" approach.^'^^ While 
the latter takes into account the lateral growth of an oblique surface, the former applies to discrete 
systems for which the nonlinear effects are negligible. Our subsequent analysis follows this 
tradition. 

I. Harmonic approach. 

The kinetic equation for the time evolution of the height within this framework is: 

^_,dh{^ = F + vV'h{x, t) - ^oy sin[7o(/i(f , t) + d{x))] + C(f , t), (2) 

where // is a microscopic kinetic coefficient; F is the driving force proportional to the difference in 
chemical potential for solidification or determined by the rate of deposition for epitaxial growth; 
u is the surface tension; y is the coefficient of the leading term due to the discreteness (higher 
harmonics are less relevant^) and Jq — ^ where b is the vertical lattice spacing; C(^)^) is the 
noise term with {C{x, t)C{x', t')) = 2D5{x — x')5{t — t'). Note that this equation of motion can be 
derived from an Hamiltonian //~^^ = + C) detailed balance is obeyed, and at (or slightly 
off) equilibrium the temperature of system is T = D/i. Eq.(2) controls the lattice growth upon 
any generic rough surface, whatever the origin of the imperfections is. We consider the finite 
quenched fluctuations d{x) of the substrate height (for a schematic description see Fig. 1.) to 
have only short range correlations and to be at least of the order of the lattice spacing of the 
bulk solid b (For \d{x)\ « b the disorder is irrelevant). Defining 9{x) — ^od{x), the phase-like 
variables 9{x) obey (^e^^i^)e^'^(^(^')^ = g^'^^'^(^x — x'), where a is the lattice spacing in the horizontal 
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planes. 

In absence of disorder {d{x) = 0), lattice effects have been studied within the harmonic 
framework by Chui and Weeks^ (CW) and later by Nozieres and Gallet'' (NG). Their important 
findings were: For D/x = T > T^, the kinetic (as the static) behavior is unaffected by the lattice 
potential C{L,t) ~ (In L)/(r/L^) (corresponding to a = = 0, z = 2) and the macroscopic 
mobility jJiM = liniiT'^o ( =< f| >- the average growth rate) is finite. For T < T^, in the 
smooth phase, the surface tends to be pinned at the periodic minima, and /lAf = (with a finite 
jump at T = Tr). For finite F, therefore, the growth is "activated" with nucleation^^'^° of higher 
" islands" . 

To apply the renormalization group (RG) approach to Eq. (0) with the substrate disorder, we 
use the Martin-Siggia-Rose (MSR) formalism.^^ In addition to h{x,t), an auxiliary field h{x,t) 
is introduced as well as their conjugate "sources" J{x, t) and J{x, t). The thermal noise and the 
quenched disorder are averaged upon to yield the following averaged generating functional: 



Z[J, J] = J VhVhexp{ J J £xdt[Jh + Jh + D^^h^ - h{— - /xz/V^/i)] 

d'^xdtdt'h{x, t)h(x, t') cos[y{h{x^ t) — h(x, t'))]} 



+ I I I d^xdtdt'Vh{x,t)Vh{x,t'), (3) 



2a' 
2' 

with g = y'^a'^. The last term has been included since it is generated under renormalization. 

Our RG analysis of this effective field theory follows closely that of Goldschmidt and Schaub 
for the XY model with random anisotropies^^ (details of the calculation will be published 
elsewhere^^ ) . Under rescaling x — > e'x and t e'^t ( Z = ln6 where h is the rescaling fac- 
tor) the RG analysis yields, to lowest nontrivial order in the following recursion relations: 



dl 

(4) 

The last equation provides the other parameter of the expansion 5 — — 1 (the deviation 
from the critical point) and these equations are of the first order in S. 7 remains at its bare 
value 7 = 7o and u at uq (we may rescale h such that = 1). We first look at the static 
properties by defining T — Dji. It obeys ^ = and the last equation can be written as: 
^ = 2(1 — Y~)g — ^g'^ with Tgr = We recognize the static equation of Cardy and Ostlund 
for the random-anisotropy XY model. Toner and DiVincenzo have analyzed these equations in 
the context of equilibrium crystal surfaces with bulk disorder in a limit where their bulk disorder 
is equivalent to the substrate disorder considered here.^^ So their results directly apply: For 
T > Tsr, — > 0, the disorder is irrelevant, and w{L) ~ (InL)^/^. For T < Tgr, g approaches 
a finite value g* ~ —5 = 1 — tJ^, disorder is relevant, and the surface becomes superrough^^ 
w{L) ~ In L. 

We now turn to study the kinetic properties. In the high-noise regime {S > 0) we find 
z = 2. In the low- noise regime z increases (the width spreading becomes slower) continuously as 
z — 2 + 4:^/c\5\ with c = |e^^ ~ 0.7931, where E is the Euler constant. The mobility hm far from 
the transition in the high Dii{= T) noisy phase assumes a finite value. Approaching Tgr from 
above, however, vanishes with 5 continuously. Integration of the recursion relation yields 



= 2^1^!'' '^ith rj = ~ 1.58. In the low-noise phase hm vanishes as |; ~ |^e""^^*~^^' ~ 
g-4V3|5|z fpjjg behavior of z and /im as function of T is summarized in Fig. 2. 

Physically the dynamic behavior in the superrough phase may be understood as follows: 
the superrough phase is a phase in which the disorder barely dominates over the thermal 
fluctuations.^^ The increased roughness is the result of the surface attempt to balance the ten- 
dency of the surface to adjust to the substrate (making h + d an integer multiple of b) without 
paying too much in elastic energy. Although the pinning is not uniform as with a smooth sub- 
strate the effect of locally preferred locations is enough to slow the spreading of the "surface 
width " (as manifested by z > 2) and to prevent it from moving with a uniform average velocity 
when an infinitesimal driving force is F applied. Naturally if a finite force is applied the sur- 
face will move, on the average, at a constant velocity. This motion will wipe out the pinning 
effect (as it does to the periodic potential in absence of disorder). If the force F is small the 
behavior described here will hold up to a scale L < L* ^ aF~^/^ and the effective mobility will 
be = InL*) ^ /xo(^*)"^^'^' ~ i^oF'^^^^^. It should be emphasized that this yields another 
important finding namely the non-linear response to a small F in presence of which the averaged 
velocity scales as F^v^l'^l^^. This explicit behavior was derived based on scaling near the critical 
point but is also consistent with activated dynamics over free-engery barriers^^ -E'(^) given by 
e{T) InL with e{T) = 4:y/^{Tsr - T). 

II. The "non-linear" approach. 

Kardar, Parisi, and Zhang (KPZ) have pointed out that when the lateral growth of an 
oblique surface is accounted for the most relevant effect is the addition of a term of the form 
|(V/i)^ to the growth equation.^ In absence of any lattice or disorder effects the KPZ equation 
is A*~^f| — F -\- i/V^/i -I- |(V/i)^ -I- C(^)^)) and the exponent a and z change^ from a — Q and 
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z = 2 for A = 0, to a ~ 0.4 and z ~ 1.6. 

The question of whether a phase transition may occur in presence of nonhnearity and a lattice 
with a perfectly flat substrate has been considered in a number of recent simulation of deposition 
(or growth) of discrete particles. ^"-"^^ E.g., the observed transition between logarithmic and power- 
law behavior of w^iV) by Amar and Family^ has been attributed to an effective vanishing of the 
nonlinear term in their discrete model. An extensive analytic study which includes both 
the lattice effects and the nonhnearity was performed by Hwa, Kardar and Paczuski^^ (HKP). 
Studying the intermediate scale L < L*, HKP identified two phases. One is the high temperature 
(strong noise) rough phase crossing over to KPZ scaling. Approaching the transition from this 
phase the mobility vanishes as (In jT — Tcl)"^ ■ They also argue that for T <Tc (lower noise) the 
surface is flat. This identification requires some caution since the generation of a term y2 cos(^/i) 
from the contraction of the terms yi sin(^/i) and |(V/i)^ was not considered. Combining both 
terms into the form \y\ sin(^/i + 0{l)) (with y"^ — yf + ^(0 — ' ^® ^^^^ 

T < Tc the flow is indeed toward \y\ — > oo but the phase shift angle is rotating hke 9{l) = col 
with I, CO ~ X/iy. Thus, this low temperature phase is not characterized simply by an increase of 
the periodic potential. Higher order terms in the recursion relations of yi and y2 and that of u 
will be required to identify with more confidence the nature of the low-temperature phase. 

Turning now to the effect of substrate disorder in presence of the KPZ nonhnearity: The RG 
analysis becomes much more complex. It turns out that a systematic expansion in the parameter 

5 — 1 requires the consideration of diagrams containing up to three non-trivial loops. 

Sophisticated techniques, based on dimensional regularization, were employed to extract their 
singular parts. The following recursion relations for x — > xe' i — > te^'' were derived (7 = 70) 



^ = 2F + 7rA (5c) 
at 

f - <£-^-^)- 

(6) 

where i^o is set to be 1, and d ~ 180.08. Combining the equations for D and /x together we find 
(r = D/.): 

The most important observation is that the nonhnearity increases the effective temperature 
T(/) and, as a result, the effective couphng g{l) becomes smaUer. The only fixed point has 

A t 

T ~ T^e^ — > oo and therefore 51 — > 0. The equation for g'(Z) can be integrated exactly: 

1 =le^(0_^e^(0 (dxe-<^\ (8) 
g{l) go 87r2 Jo ' ^ ^ 

where s(a;) = — 2)x + TQ^{e^ — 1). The behavior on long scales will have g ^ ^ and 
the effective KPZ coupling K — [note that to these orders in 5, g and A^, the flow of this 
coupling is unaffected by g ] will control the behavior. 

However, since the scale associated with the increase of the KPZ coupling is exponentially 
large Lk — ae ^0 it is hkely to be larger than L* (the scale set by F) . Therefore on scales smaller 
than L* a rough surface will be observed but it will be in a crossover regime. If Aq and/or Tq 



are small (especially with qq large), g{l) will decay to zero quite slowly. This will be observable 
in the mobility which has the scale dependence = yUoexp{ — ^^^^ !odl'g{l')} with g{l) given 
in Eq.(|^). ~ /^(^* = InL*) will be drastically reduced by the disorder effects (as they decay 
on intermediate scales) compared with its bare value /xq; the same will hold for the width w{l) 
when compared with the (7 = case. 

To summarize we have investigated how the scaling properties of growing crystalline surfaces 
are affected by disorder in the substrate. For T < T^t in the harmonic approach, the surface is su- 
perrough with anomalous dynamics. The height-height correlations are C(L, r) ~ (lnL)^/(L/r^) 
with z = 2 + 4v^(l - ^). At the same time the response becomes nonlinear: v ~ F^+i with 
C, = 2a/c(1 — T/Tsr)- In presence of the KPZ non-linearity the complex reduction in the width 
and the mobility of the rough surface have been obtained. These effects may be discernible in 
future precise measurements of solidification and epitaxial deposition processes. 



9 



Acknowledgments 

We are thankful to T. Hwa, M. Kardar, G. Grinstein, B. Schmittmann and especially to D. 
Huse for most useful discussions. We are also indebted to Y. Goldschmidt for bringing to our 
attention Ref. 22. This work was partly supported by a grant from the Corporate Research 
Laboratories of the Eastman-Kodak Company in Rochester. 



10 



References 

[1] For reviews of recent experimental and theoretical developments see e.g., Kinetics of Order- 
ing and Growth at Surfaces, edited by M. Lagally (Plenum New York, 1990), Dynamics of 
Fractal Surfaces, edited by F. Family and T. Vicsek (World Scientific, 1991) and references 
therein. 

[2] F. Family and T. Vicsek, J. Phy. A 18, L75 (1985) 

[3] S. Edwards and D. R. Wilkinson, Proc. R. Soc. Lond. A 381, 17 (1982) 

[4] M. Kardar, G. Parisi, and Y. C. Zhang, Phys. Rev. Lett. 56, 889 (1986) 

[5] L.H. Tang and T. Nattermann, Phys. Rev. Lett. 66, 2899 (1991) 

[6] S. T. Chui and J. D. Weeks, Phys. Rev. Lett. 40, 733 (1978) 

[7] P. Nozieres and F. Gallet, J. Phys. (Paris) 48, 353 (1987) 

[8] J. G. Amar and F. Family, Phys. Rev. Lett. 64, 543 (1990) 

[9] H. Yan, D. Kessler, and L. M. Sander, Phys. Rev. Lett. 64, 926 (1990) 

[10] H. Guo, B. Grossman, and M. Grant, Phys. Rev. Lett. 64, 1262 (1990) 

[11] B. M. Forrest and L. Tang, Phys. Rev. Lett. 64, 1405 (1990) 

[12] Y. P. Pellegrini and R. JuUien, Phys. Rev. Lett. 64, 1745 (1990) 

[13] T. Hwa, M. Kardar, and M. Paczuski, Phys. Rev. Lett. 66, 441 (1991) 

[14] S. T. Chui and J. D. Weeks, Phys. Rev. B 14, 4978 (1976) 

11 



[15] Substrate disorder was found to be relevant in models of sand piles by J. Toner, Phys. Rev. 
Lett. 66, 679 (1991) 

[16] J. Krug and H. Spohn, Phys. Rev. Lett. 64, 2332 (1990) 

[17] J. M. Kim, T. Ala-Nissila, and J. M. Kosterlitz, Phys. Rev. Lett. 64, 2333 (1990) 
[18] D. A. Huse, J. G. Amar, and F. Family, Phys. Rev. A 41 , 7075 (1990) 
[19] W. van Saarloos and G. H. Gilmer, Phys. Rev. B 33, 4927 (1986) 

[20] For the effect of bulk disorder on the growth rate of vicinal surfaces see J. Toner, Phys. Rev. 
Lett. 67, 2537 (1991) 

[21] P. C. Martin, E. Siggia, and H. Rose, Phys. Rev. A 8 423 (1973) 

[22] Y. Y. Goldschmidt and B. Schaub, Nucl. Phys. B 251 77 (1985) 

[23] Y. C. Tsai and Y. Shapir - in preparation 

[24] J. L. Cardy and S. Ostlund, Phys. Rev. B 25, 6899 (1982) 

[25] J. Toner and D. P. DiVincenzo, Phys. Rev. B 41, 632 (1990) 

[26] Implications to the dynamics of charge-density waves and vortex glasses will be discussed 
in Ref (23). 



12 



Figure Captions 

Fig. 1. - A two dimensional cut (along a lattice plane perpendicular to the disordered substrate) 
of the three dimensional system. 

Fig. 2. - The dependence of the hnear response macroscopic mobility /xm (bold hne) and the 
dynamic exponent z (dashed line) on temperature for the harmonic model. The arrows indicate 
the appropriate scales on the vertical axis. (R - the rough phase T > T^^, SR - the superrough 
phase T < T^r) 
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